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Abstract 

The recursive method for computing the generalized LM- inverse of a constant rectangular ma- 
trix augmented by a column vector is proposed in |16l I17| . The corresponding algorithm for the 
sequential determination of the generalized LM-inverse is established in the present paper. We 
I prove that the introduced algorithm for computing the generalized LM inverse and the algorithm 

for the computation of the weighted Moore-Penrose inverse developed by Wang in [23] are equiv- 
alent algorithms. Both of the algorithms are implemented in the present paper using the package 
MATHEMATICA. Several rational test matrices and randomly generated constant matrices are tested 
^ ■ and the CPU time is compared and discussed. 

00 
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1 Introduction 

As usual, let C be the set of complex numbers, C™^" be the set of m x n complex matrices, and 
C;"^" = {X e C™^" : rank(X) = r}. For any matrix A e C™^" and positive definite matrices M and 
N of the orders m and n respectively, consider the following equations in X , where * denotes conjugate 
, and transpose: 

(1) AXA^A (2) XAX = X 

(SM) (MAX)* = MAX (AN) {NXA)*=NXA. 

The matrix X satisfying equations (1), (2), (3A/) and (AN) is called the weighted Moore-Penrose 
inverse of A, and it is denoted hy X = A\j ^. Especially, in the case M = /,„ and N = In, the matrix 
X = A\,j ^ becomes the Moore-Penrose inverse of A, and it is denoted hy X = A^ . 

Various methods for computing the Moore-Penrose inverse of a matrix are known. The main methods 
are based on the Cayley-Hamilton theorem, the full-rank factorization and the singular value decom- 
position (see the example [1]). The Greville's partitioning method, introduced in [4], is one of the most 
efficient algorithms for computing the Moore-Penrose inverse. Two different proofs for the Greville's 
method were presented in [21 124j . Udwadia and Kalaba gave an alternative and a simple constructive 
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proof of Grcville's formula in [19]. In [3] Fan and Kalaba determined the Moorc-Pcnrose inverse of 
matrices using dynamic programming and the Belman's principle of optimality. Sivakumar in [12) used 
the Greville's formula for Al = [Ak-i |afc]^ and just verified that it satisfies the four Penrose equations. 
This provides a proof of the Greville's method by the verification. 

The Grcville's algorithm is used in various computations, where its dominance is verified over various 
direct methods for the pseudoinverse computation. The computational experience presented in [7] 
is: "When applied to a square, fully populated, non-symmetric case, with independent columns, the 
Greville's algorithm was found that the approach can be up to 8 times faster than the conventional 
approach of using the SVD; rectangular cases are shown to yield similar levels of speed increase". The 
Greville's method has been used as a benchmark for the calculation of the pseudo-inverse. 

Due to its computational dominance, this method has been extensively applied in many mathemat- 
ical areas, such as statistical inference, filtering theory, linear estimation theory, optimization and more 
recently analytical dynamics [20] (see also [6]). An application in a direct approach for computing the 
gradient of the pseudo-inverse is presented in [7]. It has also found wide applications in database and 
the neural network computation [5] . In the paper |S] , the sequential determination of the Moore-Penrose 
inverse by dynamic programming is applied to the diagnostic classification of electromyography signals. 

There is a lot of extensions of the partitioning method. Wang in [23] generalized Greville's method 
to the weighted Moorc-Pcnrose inverse. Also, the results in |23[ are proved by using a new technique. 
Udwadia and Kalaba developed the recursive relations for the different types of generalized inverses 
[2Tl[22j. Finally, the Greville's recursive principle is generalized to various subsets of outer inverses and 
extended to the set of the one- variable rational and polynomial matrices in [15j . 

The algorithm for the computation of the Moore-Penrose inverse of the one-variable polynomial 
and/or rational matrix, based on the Greville's partitioning algorithm, was introduced in |13| . The 
extension of results from |13[ to the set of the two- variable rational and polynomial matrices is introduced 
in the paper [TO] . 

The Wang's partitioning method from [53], aimed in the computation of the weighted Moore-Penrose 
inverse, is extended to the set of the one- variable rational and polynomial matrices in the paper [14| . 
Also the efficient algorithm for computing the weighted Moore-Penrose inverse, appropriate for the 
polynomial matrices where only a few polynomial coefficients are nonzero, is established in [9]. 

In the paper [6] the authors derived a formula for the computation of the Moore-Penrose inverse of 
M*M and obtained sufficient conditions for its nonnegativity, where M ~ [A\a]. 

On the other side, there are a few articles which are interested in with computation of the generalized 
LM-inverse. The definition of the LM-inverse and the recursive algorithm of the Greville's type (for a 
matrix augmented by a column vector) are given in |16 [ I17[. The recursive relations in |16[I17[ are proved 
by direct verification of the four conditions of the generalized LM-inverse. Also, these formulae are 
particularized to obtain recursive relations for the generalized L-inverse of a general matrix augmented 
by a column |17j . The recursive relations for the determination of the generalized Moore-Penrose 
Af-invcrse are derived in |18| . Separate relations for the situations when the rectangular matrix is 
augmented by a row vector and when such a matrix is augmented by a column vector arc considered 
in [18] . The alternative proof for the determination of the generalized Moorc-Pcnrose Af-invcrse of 
a matrix through the direct verification of the four properties of the Moore-Penrose Af -inverse are 
presented in [TT] . 

It is not difficult to verify that the conditions which characterize the generalized LAf-inverse are 
equivalent with the corresponding equations characterizing the weighted Moore-Penrose inverse. More- 
over, the matrix norms minimization used in (3) and (4) in the article |16j also characterizes the weighted 
Moore-Penrose inverse. Therefore, the generalized LAf-inverse and the weighted Moore-Penrose inverse 
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are identical. In the present paper we compare the corresponding algorithms. It is realistic to predict 
that algorithm for computing the weighted Moore-Penrose inverse from |23| and the algorithm for the 
computation of the generalized LM-inverse, introduced in the present paper and based on the results 
from |161 I17j are the same. Verification of this prediction is the main result of the present paper. 
Therefore, the present paper is continuation of the papers [HI HH UHl E] • 

The structure of the present paper is as follows. In the second section we restate the representation 
of the generalized LM -inverse from |16[ llTj as well as the representation and algorithm for computing 
the weighted Moore-Penrose from [53] . We also introduce an effective algorithm for construction of the 
generalized LM-inverse directly using its representation proposed in |161 117] . In the third section we 
provide a proof that two algorithms from the second section are equivalent. Implementation of both 
algorithms and a few illustrative examples are presented. 



2 Preliminaries and motivation 

The recursive determination of the weighted Moore-Penrose inverse ^]^^ ^ is established in |23| . 

Let A S C™^" and Ak be the submatrix of A consisting of its first k columns. For fc = 2, . . . , n the 
matrix A^ is partitioned as 

Ak = [Ak-i I ak] , (2.1) 

where ak is the fc-th column of A. 

Theorem 2.1 (G.R. Wang, Y.L. Chen [53]). Let A e C™^" and A^ be the submatrix of A consisting of 
its first k columns. For k ~ 2, . . . ,n the matrix A^ is partitioned as in l\2.1\i . and the matrix G C'^^'^ 
is the leading principal submatrix of N , and Nk is partitioned as 



Nk 



Let the matrices Xk-i and Xk be defined by 



Nk-i h 



nkk 



Xk-i = {Ak-i)\ 



Xk = {Ak)\ 



the vectors dk, Ck be defined by 



dk = Xk-iOk 

Ck = ak - Ak-idk = {I - Ak-iXk-i)ak. 



Then 



Xk 



Xk-i-{dk + (/ - Xk^iAk^i) N^.^lkM 
bt 



whe 



{clMck)-^clM, Ck^Q 
{dlNk-i - ID Xk^i, Cfc=0, 



and 5k = Ukk + dlNk-idk - (d^k + l*kdk{s)) -~ll{L ~ Xk-iAk-i) N^^h- 

According to the above theorem the next Algorithm 12.11 is introduced in [23] . 



(2.2) 

(2.3) 

(2.4) 
(2.5) 

(2.6) 
(2.7) 
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Algorithm 2.1 Computing tlie wciglitcd M-P inverse A\j using algorithm from f23| . 

Require: Let A € C™^", M and N be p.d. matrices of the order m and n respectively. 

1: Ai = ai. 

2: if oi = 0, then 

3: Xi = {alMaiy^alM; 

4: else 

5: Xi = 0. 

6: end if 

7: for k = 2 to n do 

8: dk=Xk-iak, 

9: Ck = ak - Ak-idk, 

10: if Cfe 7^ 0, then 

11: bl = (cp/cfc)-i4M, goto StepE^ 

12: else 

13: Sk = ukk + dlNk^idk - {dlh + l*kdk) - l*k{I - Xk^iAk-i)N-\lk, 

14: bl = S^\dlNk-i-ll)Xk-u 

15: end if 

^ _ r Xk-i - (dk + (/ - Xk-iAk-i)N^\lk)b*k ' 

16: Xk-^ 

17: end for 

18: return At = Xn- 



Also, the next auxiliary Algorithm 12.21 required in Algorithm 12. II is stated in [23] . 



Algorithm 2.2 Computing the inverse matrix N 



Require: Let A^^ 



Nk- 
II 



h 

nkk 



g C be the leading principal submatrix of p.d. matrix N . 



1: iVf 1 

2 

3 
4 
5 



6: 



for fc = 2 to n do 

9kk = {nkk - ^I■^fcl\^fc)~^ 

/fe = —9kkNf^\lk, 

Ek-i=N-\+g-^fkrk, 
Ek-i fk 

fk 9kk 



7: end for 

8: return N^^ = N^^^. 



The definition of the generalized LM- inverse is given in [16], and it is based on the usage of the 
linear equation Ax = b, where A is an m x n matrix, b is an m-vector and x is an n-vector. The matrix 
A]^j^j- is such that the vector x, uniquely given by a; = A^^j^^jb, minimizes both of the following two vector 
norms (conditions (3) and (4) from [16]) 

G=|i L'^'{Ax-b) f=\\ Ax-b\\l 
H^\\M'/'x \M\x\\l,, 

where L is an m x m symmetric positive-definite matrix and M is an n x n symmetric positive-definite 
matrix. 
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The recursive formulae for determining the generahzed LM-inverse of any given matrix A are 

introduced in |161ll7j . and they are restated here for the sake of completeness. 



Theorem 2.2 (F.E. Udwadia, P. Phohomsiri |16[|17j ). The generalized LM-inverse of any given matrix 
B ^ [A \ a] Cz C'"^" is determined using the following recursive relations: 



bIm 



[A\a 



t 



p4 



I -AA 



h 



L,M. 



a = 0, 



(2.8) 



where A is an m x [n — 1) matrix, a is a column vector of m components, 



and 



_ d^L 
^ ~ d^Ld' 

Olxm 

P = 



4- 

u = 



q^MU 
, 1 



V + p 
-1 



, V ^ A 



L,M- 



Note that L is a symmetric positive definite m x m matrix, and 

M fh 



M = 



— 7^ 

m m 



(2.9) 



where M is a symmetric positive-definite nxn matrix, M_ is a symmetric positive-definite (n—l) x (n — l) 
matrix, fh is a column vector o/ n — 1 components, and rh is a scalar. 

Theorem 12.21 assumes in i? = | a] that the matrix B is obtained augmenting the matrix A by 
an appropriate column vector a. In the rest of the paper we assume that i? = [A | a] is just the 
partitioning (|2.ip : B = Ak, A = A^-i, a = Uk- Moreover, it is clear that the following notations 
immediately follows from Algorithm 12. II 



bIm 



Xk, A^ 



Xk-i- 



Also, we use the following denotation for the matrix M defined in 

A/fe-l TOfc 



M 



ruk.k 



(2.10) 



Finally, the vector d corresponding to the first k columns of A is denoted by d^- 
According to the above Theorem 12.21 we introduce the next algorithm. 
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Algorithm 2.3 Computing tlie Lil/-inverse using the representation from 



Require: Let A e 



L and M be p.d. matrices of order m and n respectively. 



1 
2 

3 

4 
5 
6 
7 
8 
9 
10 
11: 
12: 

13: 

14: 

15: 
16: 

17: 



Ai = ai. 

if Ci = then 

Xi = {alLai)-^alL 
else 

Xi = 0. 
end if 

for fc = 2 to n do 

dk = {I — Ak-iXk-i) ak, 
p ^ [I - Xk-iAk-i) M-} 
i^ dk^O then 

hi = goto Stepm 



rrik, 



-'k 
else 



q = 

u = 
K = 

end if 

Xk - 



Xk-iak 
-1 

Xk-i ' 




p 



q*AIkq 
Xk-l 



18: end for 
19: return A\ 



L.M 



MkU, 

Xk-ittkbl 

K 

Xn. 



It is clear from restated definitions that the generalized _Li\f -inverse is just the weighted Moorc- 
Penrose inverse. Therefore, Algorithm 12.31 and algorithms 12.11 12.21 together produce identical result - 
the weighted Moore-Penrose inverse of the given m x n matrix. In the next section we compare the 
described algorithms. 



3 Comparison of algorithms 

Theorem 3.1. Aloorithm \2.3\ is equivalent to aloorithms \2.1\ and \2.SX 

Proof. In order to ensure unambiguous, during the proof we assume that the symbol in a superscript 
denotes terms from the Wang's algorithm; similarly we use the convention that U , as a superscript, 
denotes terms from the Udwadia's algorithm, elsewhere it is necessary. Since the iM-inverse is just the 
weighted Moore-Penrose inverse, we conclude that the matrix M in Algorithm 12. II is just the matrix L 
in Algorithm 12.31 and the matrix N in Algorithm 12. H is analogous with the matrix M in Algorithm 12. 31 
Therefore, it is not necessarily to mark the matrices L, M, N and A by appropriate superscript. 

We prove the theorem by verifying the equivalence of the outputs from the corresponding algorithmic 
steps of mentioned algorithms. The proof proceeds by the mathematical induction. 

The proof for the case fc = 1 in view of Step 3 in both algorithms is trivial. Assume that the 
statement is valid for the first k — 1 columns, i.e. 

Xt, = X^_, = = 4Af.- (3.11) 
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Now wc verify the inductive step. Wang used the matrix Xk in the form 

xr-i-{d^ +(1- x^UAk-i) N^\h){bir 



while Udwadia observed two cases, as in 
If we denote with 



then the equahties in (|2.8p become 



(3.12) 



(3.13) 



(3.14) 



Let us show that the output of Step 9 from Algorithm 12.11 is the same as the output of Step 8 from 
Algorithm [ 



Cfc 



W — 



Now we show that Step 11 from Algorithm 12.11 and Step 11 from Algorithm (53] are equivalent. As it is 
stated above = d,^ > so that in the case 7^ we have 



*\U 



In a similar way it can be verified that Step 14 from Algorithm 12.11 is equivalent to Step 15 from 



Algorithm 12.31 In the case = we can start from the statement in Step 15 from Algorithm 

{blf = {q*MuU)l{q*Muq). 
From Step 9 of Algorithm 12.31 and the inductive hypothesis the following holds 
p = (/ - X}!_^Au-i) N^\h = {I- X]^_,Ak-i) N^\h, 
so that we derive the following: 



[iXj!_,ak+pr I -1] 



Nk-i Ik 



nkk 



XU 




^k-i 




= [{Xj!_,ak+prNk-i~ll I {X^_,ak+prik-nkk] 

= {X^_,ak +prNk^iXj!., - IIX^_, {since Xj^_,ak = X^^au = df} 

= {dirN,^^Xr^,-llX^^,+p*Nk-^Xr^, 

= {{dirNk-^ II) Xj^., + ((/ - Xr.,A,.,) N^l,l,y N,^,X^_, 

= mr^k-, - II) x^_, + [iiN-\ (/ - xjr.,A,.,)*) Nk-^xz, 

= {{dirNk-^ II) X^_, + IIXY., - llN^l, {X^.^Ak-iY Nk^.XY^,. 

Since A^fc_i is the symmetric positive definite applying equality (4iV) together with p. lip the following 
holds 



{Xl^_^Ak-iy Nk-i = (Nk-iX^^^Ak-i) " = Nk-iX^_^Ak-u 
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and later 





-11) 




7* vW 




= {{dirNk-, 


-11 






- x^:,A,^,xz,) 


= {{dirNk-i 


-11 








= {KrNk-i 


-11 









Moreover, we have 



q'^Mkq = [iXj^_,ak+pr \ - l] 



It 



Ik 

nkk 



Xj^_,ak+P 
-1 



= [{XlUak + p)* Nk-i - l*k I {Xj!_^ak+pyik-nkk] 



Xj!^,ak+P 
-1 



iiXj^_,ak+prNk-i - ll){X^_^ak+p) - {X}^_^au +prik + n^k 
{X]i_^ak+prNk^^{X^^^ak + p) ~ ll{Xj!_^ak + p) - {Xj!_,ak +prh + n^u 



(since {X\!^^ak = X^_^ak = < }) 
{df +pTNu-r{df +p) mf+p) {dl 



■w 



-pTh 



' nkk 



(dl)'^ Nk-id"^ 



■■k) JVfc-iu;, - lldY - [dlY^h + nkk - l*kP 
+ p*Nk-idY + {dD'^Nk-iP + p*Nk^ip - p*lk. 



Furthermore 



p*Nk-id^ + idl)'^Nk-ip + p*Nk-ip-p*lk = 0. 



(3.15) 



First we show that p*Nk-idj: = 0, as follows 



p*Nk-idY = {{I - Xl''_,Ak-^) N-},lkY Nk^iX^^.ak 

= llN-\ {I - X]^_,Ak-iY Nk-iXY_,ak 

= IIN^\ (Xk-i - [xf^^Ak-iY Nk-i) X^_,ak 

= IIN~\ (iV.^i - Nk^iXY^^Ak-i) XZiak 

= llN^\Nk-i {I - xZiAk-i) X^^.ak 

= ll{l-Xr^,Ak-i)xr^,ak 

= I'k {X^_;i^ - X^_-i^Ak-lX^_l) ttk 

= ll{X]^:,-Xr.,)ak^O. 



Also, from the above equality, we have 

{dirNk-ip = {p*Nk-idYy 



Finally, the last term of the sum in the left hand side of p.lSp . is equal to p*Nk-ip — p*lk, and it is 



About the generalized LM-inverse and the Weighted Moore-Penrose inverse 



9 



also equal to zero: 

p*Nk-ip - p*lk = P* {Nk-i {I - X]^^,Ak-i) N^\lk - Ik) 
- p*{-Nk-iX^_,Ak-iN^},lk) 
= p* {-ixZ,Ak-ir Ik) 

= ((/ - N^\lk)* {-AU,iX*,^,rk) 

= -llN^\Al_,{Xl_,)'^h + llN^l,Al„,{XLi)'^AU{XUrh 

= ~llN-\Al_,iXl_,rh + llN-},AU,iXUrh - 0. 

Continuing the transformation for Alj-q we have 

q^Mkq = {dirXk^^d^ - lld^ - {dl^h + nuk X^iA,.,) N'^h 

= 4. 

Now we are able to continue the rest of our proof. According to the equality 

(6:.)^ = {q*AhU)/{q*A'hq), 
Step 15 from Algorithm 12.31 produces the result 

{bif = s^\diN,^,-i:)xz„ 

which is identical to the output (6^)^, derived in Step 14 from Algorithm 12. II 



According to Step 16 of Algorithm 12.11 and Step 17 of Algorithm 12.31 the generalized LM-inverse 
and the weighted Moore-Penrose inverse of the first k columns of A are identical: 



Xt = 



Xj!_,-{X^_,a,+p) [hir 

'l-XZ,Au^i)N^\lk) 
If 



vW 
^k-1 



ibl)''' 



(3.16) 



Xl 



Finally, for the case fc = n, it immediately follows that A^j^ 
from both algorithms are identical. □ 



A\ r which means that the outputs 



4 Examples 

In order to compare the algorithms from the second section it is necessary to use the precise imple- 
mentation of the corresponding algorithms. Details concerning the implementation of the partitioning 
algorithm corresponding to the weighted Moore-Penrose inverse can be found in |14j . In order to com- 
pare the mentioned algorithms we developed a MATHEMATICA code for the implementation of Algorithm 
12.31 We later tested results on different types of matrices. Since the language MATHEMATICA admits 
symbolic manipulation with data, developed implementations are immediately applicable to the rational 
and polynomial matrices. 
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Example 4.1. Consider the test matrix Anxio from [25], in the ease a = 1 



- 1 


2 


3 


4 


1 


1 


3 


4 


6 


2 - 


1 


3 


4 


6 


2 


2 


3 


4 


5 


3 


2 


3 


4 


5 


3 


3 


4 


5 


6 


4 


3 


4 


5 


6 


4 


4 


5 


6 


7 


6 


4 


5 


6 


7 


6 


6 


6 


7 


7 


8 


6 


6 


7 


7 


8 


1 


2 


3 


4 


1 


3 


4 


1 


1 


3 


4 


6 


2 


1 


2 


4 


6 


2 


2 


3 


4 


5 


3 


1 


3 


4 


5 


3 


3 


4 


5 


6 


4 


2 


3 


5 


6 


4 


4 


5 


6 


7 


6 


3 


4 


. 6 


7 


6 


6 


6 


7 


7 


8 


4 


5 . 



and randomly generated symmetric positive definite matrices Lioxio and Muxn: 



L = 



280 


-5 


-133 


-27 


-12 


-93 


-133 


-42 


84 


52 


-5 


216 


-93 


-23 


141 


-2 


-108 


-48 


21 


165 


-133 


-93 


336 


1 








9 


— 5 


91 


81 


-27 


-23 


1 


260 


-62 


-42 


-40 


6 


85 


-8 


-12 


141 





-62 


278 


-63 


-19 


-135 


9 


99 


-93 


-2 





-42 


-63 


238 


68 


-27 


-80 


-34 


-133 


-108 


9 


-40 


-19 


68 


290 


12 


-233 


-244 


-42 


-48 


-5 


6 


-135 


-27 


12 


209 


-15 


-87 


84 


21 


91 


85 


9 


-80 


-233 


-15 


332 


145 


52 


165 


81 


-8 


99 


-34 


-244 


-87 


145 


402 



452 


91 


-186 


-97 


-161 


68 


28 


16 


151 


41 


-65 


91 


413 


-74 


-119 


-317 


-41 


-12 


67 


180 


-53 


54 


-186 


-74 


497 


-136 


78 


-208 


-175 


-120 


9 


-99 


6 


-97 


-119 


-136 


371 


157 


154 


129 


29 


-102 


-16 


-96 


-161 


-317 


78 


157 


444 


28 


-39 


-201 


-165 


3 


43 


68 


-41 


-208 


154 


28 


509 


52 


55 


90 


179 


-81 


28 


-12 


-175 


129 


-39 


52 


454 


-38 


-157 


145 


22 


16 


67 


-120 


29 


-201 


55 


-38 


408 


-9 


-16 


-100 


151 


180 


9 


-102 


-165 


90 


-157 


-9 


257 


-32 


14 


41 


-53 


-99 


-16 


3 


179 


145 


-16 


-32 


376 


-15 


-65 


54 


6 


-96 


43 


-81 


22 


-100 


14 


-15 


339 



The generalized LM -inverse A]^ j^j from [1^1 [H] and the weighted Moore- Penrose inverse A\j from 
[23] are both equal to 



0.755 


-0.156 


-1.917 


0.823 


0.143 


0.033 


0.383 


-0.33 


0.213 


-0.802 


0.67 


-0.542 


-0.078 


1.683 


-0.544 


-0.27 


0.003 


-0.432 


0.673 


0.075 


0.087 


-0.347 


0.346 


-0.194 


-0.881 


0.751 


-0.187 


0.002 


0.32 


-0.149 


0.711 


-1.749 


1.003 


-0.049 


0.558 


-0.724 


0.145 


0.065 


0.007 


0.128 


-0.258 


-0.245 


0.481 


-0.141 


-0.454 


-0.013 


1.181 


-0.79 


0.188 


0.109 


-0.3 


0.057 


-0.563 


1.498 


-0.907 


-1.188 


-0.468 


2.87 


-0.203 


-0.792 


-0.107 


0.233 


0.19 


1.258 


-2.561 


1.083 


0.701 


0.364 


-2.009 


0.492 


0.282 


0.009 


0.395 


-0.407 


-0.481 


0.786 


-0.192 


0.472 


0.156 


-0.533 


-0.675 


0.507 


-0.02 


-0.62 


-0.008 


-0.822 


2.174 


-0.848 


-0.415 


-0.481 


1.824 


-0.22 


-0.377 


0.002 


0.013 


0.247 


0.615 


-1.143 


0.257 


0.328 


0.156 


-1.388 


0.469 


0.362 


-0.02 


0.028 


-0.008 


-0.47 


0.526 


-0.2 



The Moore- Penrose inverse can be generated in the case L = M = al, M ^ N ^ (51 [TB], and it is 
equal to 



0.294 


-0.169 


-1.511 


1.415 


-0.391 


0.041 


-0.04 


-0.26 


0.454 


-0.264 


0.23 - 


-0.067 


-0.074 


1.227 


-1.064 


0.23 


0.001 


-0.192 


0.649 


-0.103 


-0.191 


-0.102 


0.227 


-0.179 


-0.692 


0.705 


-0.214 


-0.009 


-0.254 


-0.238 


1.514 


-1.319 


0.449 


-0.165 


0.547 


-0.657 


0.376 


-0.116 


0.015 


0.179 


-0.196 


-0.456 


0.531 


-0.104 


-0.297 


-0.01 


1.035 


-0.972 


0.359 


0.108 


0.238 


0.044 


-1.065 


0.938 


-0.366 


-0.008 


-0.474 


1.663 


-1.319 


0.352 


-0.103 


-0.428 


0.224 


1.83 


-1.844 


0.414 


0.062 


0.368 


-1.349 


1.081 


-0.33 


0.006 


0.426 


-0.434 


-0.442 


0.712 


-0.156 


-0.081 


0.158 


0.029 


-0.144 


-0.034 


-0.021 


0.087 


-0.019 


-1.499 


1.448 


-0.138 


0.195 


-0.484 


1.198 


-0.791 


0.211 


0.005 


-0.19 


0.268 


0.765 


-0.906 


0.049 


. -0.081 


0.158 


-0.971 


0.856 


-0.034 


-0.021 


0.087 


-0.019 


-0.499 


0.448 


-0.138 . 
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Example 4.2. Consider the one variable test matrix 



A = 
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proposed in [25] and L (resp. M) and M (resp. N) as the identity matrices of the appropriate 
dimensions. Both of the considered algorithms produce the following Moore-Penrose inverse: 
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Example 4.3. The CPU time needed for the computation of the generalized LAf-inverse and the 
weighted Moorc-Pcnrosc inverse (according to Algorithm 12.31 and the algorithms 12.11 12.21 respectively) 
is compared in the next table. The testing is done on the local machine with the following performances: 
Windows edition: Windows Home Edition; Processor: Intel(R) Celeron(R) M CPU @ 1.6GHz; Memory 
(RAM): 512 MB; System type: 32-bit Operating System; Software: MATHEMATICA 5.2. Also, in the table 
there are arranged results obtained on the set of randomly generated test matrices Amxn and randomly 
generated symmetric positive definite matrices Lmxm (resp. M„ixm) and Mnxn (resp. Nnxn)' 



m X n 


degree 


Alaorithm\2.1\ 


Algorithra 12.31 


5x6 


1 


2.265 Seconds 


2.235 Seconds 


5x6 


2 


4.078 


4.063 


6x4 


5 


9.969 


9.625 


6x4 


10 


23.328 


21.11 


10x11 


1 


104.109 


103.484 


10x11 


2 


192.734 


186.297 


11x10 


1 


133.469 


130.125 


11x10 


2 


261.359 


227.047 



Table 1. The comparison in the efficiency on the set of randomly generated test matrices 

According to the above Table 1 results it is evident that Algorithm 12.31 produces negligibly better 
performances with respect to Algorithm 12 . 1 1 for all test cases. This fact is in accordance with the verified 
equivalence between the algorithms. 
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5 Conclusions 

Our primary idea is to show that the computational method for the generahzed LM- inverse from |16[|17| 
and the algorithms for the computation of the weighted Moore-Penrose inverse from [52] are equiva- 
lent. The effective algorithm for the computation of the generalized LM-inverse is introduced here. 
Equivalence of the considered algorithms is proved in the third section by verifying the equivalence of 
outputs generated by corresponding algorithmic steps. This paper not only compares the corresponding 
algorithms but also compares the performance of two approaches of finding the Moore-Penrosc inverse. 
In order to compare the efficiency of corresponding algorithms we developed their implementations in 
the programming language MATHEMATICA. 



6 Appendix 

Several auxiliary procedures implemented in MATHEMATICA are described at the beginning. 

TakeCol[A_, k_] := Transpose [ Take [ Transpose [A] , {k}]]; 
TakeCols[A_, k_] := Transpose [ Take [ Transpose [A] , k]]; 
TakeElement [A_, i_ , j_] := TakeCol [ Take [ A, {i}], jj; 
TakeElements [A_, i_ , j_] := TakeCols [ Take [A, i], j]; 
T[ A_] := Transpose [A] ; 

RandomPoly [ n_ , probl_, prob2_ , var_] := Module[{S, i, 1}, 
If [Random [ Real , {0, 1}] > probl , Return [0]]; 
S = 0; 

Do[ If [Random [Real , {0, 1}] < prob2 , 1 = 1, 1 = 0]; 

S = S + var " i *Random [Integer, {-10, 10}]*1, {i, 0, n}]; 
If [(S = 0) && (probl >= 1), S = S + 1]; 
Return [S ] ] ; 

RandomMatrix [m_ , n_ , deg_ , probl. , prob2_ , var_] := Module [{A, i, j}, 
A = TablefO, {i , 1 ,m} , {j ,l,n}]; 

Do [Do [A [ [ i , j ] ] = RandomPoly [ deg , probl , prob2 , var ] , {1,1 ,m} ] , { j , 1 , n } ] ; 
Return [A] ] ; 

Implementation of Algorithm 12.31 is given by the following function. 

AL1V!W[A., M., N.] : = 

Module [{I, m, n, d, c, tmp , tmpl , u, q, p, a, Al , 1, X, Nl , N2 , nil}, 
a = TakeCol [A, 1] ; {m, n} = Dimensions [A] ; 
If [ Together [a] = 0*a, 
X = T[a] , 

X= Inverse [T[a] .M. a] .T[a] .M ]; 
Al = {}; 
Do[I = IdcntityMatrix [ i - 1]; 

a = TakeCol [A, i]; Al = TakeCols [A, i - 1]; 

nil = TakeElement [N, i, i]; Nl = TakeElements [N, i - 1, i - 1]; 
N2 = TakeElements [N, i, i]; 1 = T [ TakeCols [ Take [N, {i}], i - 1]]; 
d = ( IdentityMatrix [m] - Al.X).a // Together; 
p=(I-X.Al).Inverse[Nl].l // Together ; 
If [Together [d] =!= 0*d, 

b = Inverse [T[d] .M.d] .T[d] .M; b=T[b] // Together, 

tmp = (X. a + p) // Together; q = Join[tmp, {{—1}}]; 

u = Append [X, Table [0, {j,m}]]; 

b = (T[q] .N2.u) Inverse [T[q] .N2.q] [ [1] ] ; 

b = T[b] // Together ; 

]; 

tmpl = X - (X.a - p).T[b] // Together; 
X= Together [ Join [tmpl , T[b]]]; 
, {i, 2, n}]; 
Return [X // MatrixForm ] ] ; 
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Implementation of Algorithm I2.1l is obtained by slightly adopting the MATHEMATICA code described 
in [H]. 

AWang[A., M., N.] : = 
Module [{I, m, n, d , c, a, Al , 1, del, X, tmp , tmpl , Nl , nil}, 
a = TakeCol [A, 1] ; 
{m, n} = Dimensions [A] ; 
If [ Together [a] ^= 0*a, 
X = T[a] , 

X = Inverse [T[a] .M. a] .T[a] .M ] ; Al = {}; 
Do[ I = IdentityMatrix [ i - 1]; 
a = TakeCol [A, i]; Al = TakeCols [A, i - 1]; 
d = X.a; c = Together [a - Al . d ] ; 
nil = TakeElement [N, i, i]; 
Nl = TakeElements [N, i - 1, i - 1]; 
1 = T[TakeCols [Take [N, {i}], i - 1]]; 
tmp = T[d] . Nl - T[ 1 ] ; 

tmpl = (I - X.Al) . Inverse [Nl] . 1 // Together; 
If [Together [c] =!= 0*c , 

b = Inverse [T[c ] .M. c ] .T[c ] .M; 

b = T[b] // Together , 

del = nil + T[d].Nl.d- {T[d].l + T[l].d) -T[l]. tmpl // Together ; 
b = T[ Inverse [ del ]. tmp .X] // Together; 

]; 

X = Together [ Join [X - (d + tmpl ).T[b], T[b]]]; 
, {i, 2, n}]; 
Return [X] ] ; 
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